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Abstract 

The Causal Dynamical Triangulation (CDT) approach to quantum gravity is a lattice ap- 
proximation to the gravitational path integral. Developed by Ambj0rn, Jurkiewicz and Loll, 
it has yielded some important results, notably the emergence of classical spacetime and short 
scale dimensional reduction. However, virtually all the results reported so far have been based 
on a single computer code. In this paper we present the first completely independent verifica- 
tion of the CDT algorithm, and report the successful reproduction of the emergence of classical 
spacetime and smooth reduction in the spectral dimension of the 2-1-1 and 3-1-1 dimensional 
spacetimes. 

1 Introduction 

The Causal Dynamical Triangulations (CDT) approach to quantum gravity has emerged as a breath 
of fresh air in the effort to quantize gravity. In CDT, developed by Ambj0rn, Jurkiewicz and Loll 
[H [21 Hj , we have an attempt at quantizing gravity that relies solely on ideas and techniques that 
have been known to physicists for quite some time — path integrals, causality, simplicial manifolds, 
finite-size scaling effects, etc. The formulation of CDT, while being extremely straightforward, has 
yielded some important results, such as the emergence of classical spacetime and dimensional 
reduction of spacetimes at short distances. 

So far, however, all the results that have been reported for CDT have either come from the 
original group headed by Ambj0rn, Jurkiewicz and Loll or from groups that have used the code or 
data sets from the original group. In this paper, we present the first independent verification of 
the CDT algorithm, using code developed from scratch in a completely independent fashion. We 
have been able to reproduce the important results reported for CDT and set the stage for further 
results, which will be reported elsewhere [HEIIT]. 

The outline of this paper is as follows. In section [2] we present a brief overview of the Causal 
Dynamical Triangulation approach. Readers interested in greater detail should refer to [2J and [3]. 
In section [3] we describe our CDT implementation, at a fairly high level. Finally, in section [4j we 
present the results from our simulation. 
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2 Causal Dynamical Triangulations 



CDT is based on the gravitational path integral 

G(gi,gf;t„t/) = I Pb]e'^^«t^] (1) 

where 

SEH[g] = I rx^{R - 2 A) (2) 

is the Einstein-Hilbert action. Here gi and gf are the spatial metrics on the initial and final slice 
and g represents all possible (spacetime) metrics that satisfy the initial and final conditions. By all 
possible metrics, we mean metrics that are not related to each other by diffeomorphisms. To define 
the action in the path integral, CDT uses Regge calculus [8J as its starting point. The basic idea 
of Regge calculus is to approximate a smooth manifold (for example spacetime) with a piecewise 
linear manifold, with curvature restricted to subspaces of co-dimension two. 

The geometrical objects used to approximate d- dimensional manifolds are d-simplices, which 
consist of {d— l)-dimensional faces, (d — 2)-dimensional hinges or "bones" (on which the curvature 
of the manifold resides), and so on down to 0-dimensional simplices (points). A d-simplex has d+1 
points. 

For manifolds without a boundary, a discrete version of the Einstein-Hilbert action can be 
formulated using [8] 

J d'^xV^R ^ 2 Yl ^h6h (3) 

hinges h 

d^x^g ^ Y.^d (4) 



to get 



SvrG ^ SttG 

h d 



where is the volume of a hinge, is the volume of a top level d-simplex, and 5h is the deficit 
angle at h. 

CDT assumes globally hyperbolic spacetimes. The existence of Cauchy hypersurfaces enables 
us to use d-simplices that span adjacent spatial slices as the building blocks for our d-dimensional 
spacetimes. A time-like {p, q) d-simplex, where p+q = d+l, has p points on the lower time slice and 
q points the adjacent higher time slice. For a space-like simplex, all points are on the same spatial 
slice. No space-like d-simplices are present in d-dimensional spacetimes. But both space-like and 
time- like types of lower dimensional simplices are present in the spacetime. The simplicial building 
blocks in CDT have an inherent causal structure built into them. To preserve this causal structure, 
CDT does not allow topology changes of the spatial slices. 

The lengths of the space-like and time- like links (1-simplices) are defined as 

^Ipace = a'^Jtime = « > (6) 

where a is the lattice spacing and a is the asymmetry parameter. In 2+1 dimensions, we have three 
types of 3-simplices, denoted by (1,3), (2,2) and (3, 1). In 3+1 dimensions, we have four types of 
4-simplices, denoted by (1, 4), (2, 3), (3, 2) and (4, 1). The volume of each type of simplex and the 
dihedral angles subtended between the faces of these simplices are given in [2]. Using these, we can 
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formulate the Lorentzian version of the Regge action in d = 3, 4 as 



space— like d—simplices 
hinges, h at h 



time— like d—simplices 
hinges,h at h 

time— like 
d— simplices ,d 

where G is the dihedral angle subtended at the hinge h. Equation ([T]) is a reformulation of ([s]) with 
space-like and time-like simplices distinguished. 
In CDT, Wick rotation is defined as [2] 

1 7 

Q 1-^ -a, a > -(d = 3), a > — (d = 4) (8) 

This leads to a pure imaginary value for the action. The gravitational path integral ([T]) in CDT is 
replaced by a sum over distinct triangulations 

ngv'^Y^^f^^^ (9) 

where C(T) is the order of the automorphism group of the triangulation. 

The Wick rotated actior|^for a triangulation in 2+1 dimensions is given by 

5(3) = i{-KoNo + K3N3) (10) 

where 



and A'o, N3 are the number of 0-simplices and 3-simplices, respectively. Equation ( 10 ) can be derived 
from Q by substituting the appropriate expressions for the simplex volumes and dihedral angles, 
using topological identities for (2+l)-dimensional Lorentzian triangulations [2j and then performing 
a Wick rotation ([s]). The topological identities eliminate all but two of the bulk variables A'^o • • • Nd- 
Substituting the Wick rotated action in ^ gives us the partition function: 

Z = YexpiKoNo-K3N3) (12) 

T 

A similar procedure in 3+1 dimensions leads to the Wick rotated action 

= i(^-(Ko + 6A)iVo + K4iV4 + A(2Aff + ivf ''))^ (13) 
and the corresponding partition function 

Z = ^exp( (ko + 6A)7Vo - - A(2A^f '^^ + ivf '^^) ) (14) 

T 



^Equations ( 10 1 and are derived for a = —1. For a ^ —1 the action will have terms involving A'^o, and 
N>/''', where A^;^^'"' and N>/'^' are the number of (3,1) and (2,2) 3-simplices. However, in 2+1 dimensions, there are 
only two independent bulk degrees of freedom, so the action still takes the form ( 10 1, although with (111 altered. See 
|2] for details. 
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where 

V3 AVE Vsfb 1,,,,, \ 

= 8G' '^^ = 768^ + 8G (2. '^'^ ^'/'^ " ') ^''^ 

and A is a complicated function that depends on a, with a = 1 corresponding to A = 0. N^'^'^^ is 

the number of (4, 1) plus (1,4) simplices, and A'^ is the number of (3,2) plus (2,3) simplices. 

To generate the distinct triangulations required for evaluating the partition function and the 
expectation value of any observable, the CDT approach uses a set of ergodic moves. An ergodic 
move applied to a triangulation yields another distinct triangulation, one that is not related to the 
first one by a diffeomorphism. Ergodicity ensures that we can go from one triangulation to any 
other triangulation by repeated application of moves. Ergodic moves also preserve the topology of 
the manifold, as well as the simplicial manifold property. The moves are described in detail in [2]. 



3 Numerical Implementation 

In our numerical implementation, each d-simplex is identified by a unique ID and is stored as a list 
(containing sublists) of points 

(ty lo hi (pO pi . . . pd) (nO nl . . . nd)) 

Here ty identifies the type of the simplex ty = 1 , 2 , . . . , d. lo and hi are the lower and higher 
time slices spanned by the simplex. The pis are the points that make up the simplex, and the 
nis are the IDs of the neighboring simplices. The simplex nk is connected on the (pO . . . pd) - 
(pk) face. The points are used only as unique labels to identify the vertices of a simplex. While we 
distinguish the points on the lower time slice from the points on the higher time slice, the ordering 
of the points in any given slice is immaterial. 

With this choice for the simplex data structure, it is fairly straightforward to implement the 
ergodic moves described in detail in [2]. It is essential that the moves maintain the simplical 
property of the manifold. This requires that two simplices should be connected across a single face 
and that a (d — l)-simplex should be shared between exactly two d-simplices. 

To generate the distinct triangulations using the ergodic moves, we apply the Metropolis algo- 
rithm [9]. The Metropolis algorithm uses a Markov chain that, upon reaching equilibrium, satisfies 
the following detailed balance condition 

P{Ti)W{Ti,T2) = P{T2)W{T2, Ti) (16) 

where P{Ti) is the probability of being in a particular triangulated geometry Tj and W{Ti,Tj) is 
the transition probability of going from triangulation Tj to triangulation Tj. The probability of 
being in a particular triangulation is proportional to the Boltzmann weig ht e-^-B. The detailed 
balance condition in this case can be written as: 

The simplest choice for the transition probability W that satisfies this condition is 

Thus a move that increases the action may be accepted or rejected, while moves that decrease the 
action or leave the action unchanged are always accepted. 
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We initialize the simulation by specifying the number of time slices T and the total volume 
of the spacetime. The initialization process starts by triangulating each spatial slice as a (d — 1)- 
sphere. The t = and the t = T slices are identified to enforce the periodic boundary conditions 
in time. This is done for technical convenience, rather than due to any physical argument. Next, 
the triangulations on adjacent slices are connected. We then apply the volume increasing subset of 
the ergodic moves until the total volume reaches the desired target volume. 

Once the simulation has been initialized, we have to run simulation for about 100,000 sweeps to 
thermalize the spacetimes. A sweep is defined as attempted moves, where A^^ is number of d- 
simplices in the spacetime. Thermalization ensures that the system is sufficiently randomized, and 
has moved away from the highly structured initial conditions. Once thermalization is completed, 
we start saving spacetimes every 100 sweeps, to build our ensembles. 

4 Results 



4.1 Phase Structure 

In 2+1 dimensions, the phase structure of the CDT model can be explored by tracking the depen- 
dence of an order parameter on coupling constants. The order parameter selected is the ratio of 
the number of (2, 2) 3-simplices to the total number of 3-simplices in the spacetime: 



N22 



N22 



Ni3 + N22 + N31 

The dependence of r on the coupling constant is shown in figure [l} 



(19) 



0.35 - 
0.3 - 



Figure 1: versus kq 

Figure [1] shows a phase transition, at Kq w 3.3. In pj, this phase transition was reported at a 
value of Kq ~ 6.6. The two results differ by a factor of two, for reasons we do not understand, but 
the qualitative results are identical. For kq > Kq, we observe that the spatial volumes of successive 
slices is not strongly correlated — spacetime has effectively decoupled into disconnected spatial 
slices. This is shown in figure [2j where we have a taken a snapshot of a simulation in this decoupled 
phase. 

A snapshot of the simulation in the other phase, characterized by kq < Kq, clearly shows the 
emergence of a well defined, extended geometry phase. This is shown in figure [3} The simulation 
fluctuates around this extended geometry, a geometry that resembles a classical universe. This clas- 
sical spacetime emerges from the simulation, and was not explicitly specified by us. In this phase. 
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Figure 2: A snapshot of a (2+l)-dimensional spacetime in the decoupled phase. 




Figure 3: A snapshot of a (2+l)-dimensional spacetime in the extended phase. 

the volume of successive spatial slices is strongly correlated, and the path integral is dominated by 
metrics that are approximate solutions of the Einstein equations. 

A thick slice is defined as the portion of spacetime that is bounded by adjacent spatial slices. 
Thick slices are labelled by half-integer values where, for example, a thick slice labelled by t = 1/2, 
is bounded by the spatial slices t = and t = 1. The volume of a thick slice is defined as the 
number of d-simplices in that thick slice. In figures [2] and [3] time is along the horizontal axis. These 
figures were generated by considering the volumes of the thick slices for half-integer values ranging 
from t = l/2tot = T — 1/2 where T is total number of time slices. We draw a circle with radius 
equal to the volume of the slice, and join circles corresponding to adjacent thick slices by linear 
interpolation. The axial symmetry in these figures is therefore a visualization artefact. Figure [4] 
shows a typical spatial slice in 2+1 dimensions, embedded in a three dimensional space [6]. 




Figure 4: A spatial slice in 2+1 dimensions. 

In 3+1 dimensions, while we are not aware of an order parameter analogous to the (2 + 1)- 
dimensional case, we do observe three distinct phases, illustrated in figures [5} [6] and [7| All three 
figures are for spacetimes with T = 64 and A'4 = 81920 

Phase A is characterized by a splitting up of the spacetime in small pieces of irregular volume, 
with each piece not extending beyond a few time slices. In Phase B, all of the spacetime volume 
collapses into a single piece of minimal extension. Both of these non-physical phases are character- 
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Figure 5: A (3+l)-dimensional spacetime in phase A, = 5.0, A = 0.0, K4 = 1.01 



ized by a small value for the ratio 



(3,2) 



A^4 



Phase C is the extended geometric phase, where there is 




Figure 6: A (3+l)-dimensional spacetime in phase B, kq = 1-6, A = 0.0, K4 = 0.65 
a strong correlation between the volumes of successive spatial slices. 



4.2 Spectral Dimension 

One of the major contributions of the CDT approach has been the introduction of the spectral 
dimension as probe of the properties of spacetime, both on short and long scales. 

The diffusion equation for a d-dimensional manifold with a smooth metric g^i, is given by \10\ [TT] 



d_ 



Kg{i,i^;a)=AgKg{i,io-,a) 



(20) 



where a is the diffusion time, A^ = —g^^V u is the Laplace operator corresponding to g^uiO and 
Kg{^,£,o] c) is the probability density of diffusion from ^0 to ^ in a time a. The diffusion processes 
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Figure 7: A (3+l)-dimensional spacetime in phase C,ko = 2.4, A = 0.6, K4 = 0.079 
we care about are ones that are initially peaked at a point ^0 



Instead of probability density, we measure the return probability, which is defined as 

^ Jm 



(21) 



(22) 



For a flat space with infinite volume, the solution of equation (20) is \10\ lllj 



e 9 



(47ra)'^/2 



(23) 



where dg{(,,^o) is the geodesic distance between ^ and ^o- The return probability for the flat space 
case is 



1 



a- 



d/2 



Taking the logarithmic derivative of the equation ( 24), we get 



a In a 



(24) 



(25) 



This definition of the spectral dimension can be extended to curved and/or finite volume spacetimes, 
with the addition of finite-size corrections, and allows one to define a generalized dimension on a 
wide variety of spaces that are not smooth manifolds [12] 

In CDT, the spectral dimension [TUl HJ [T3] of the spacetime is measured by initiating a diffusion 
process from a randomly selected simplex. A single step in a random walk is generated by moving 
to randomly selected neighbor. By initiating a large number of walks with a steps, we can calculate 
the probability of returning to the starting simplex. The return probability can be then used to 
compute the spectral dimension. 

The spectral dimension vs step size plot for a (2+l)-dimensional spacetime is shown in figure 
[8) We observe that the spectral dimension demonstrates a smooth transition from a value of about 
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3.2 h 




Figure 8: Spectral dimension for a (2+1) spacetime vs a, N3 = 80000, kq = 1.0 

2.4 at short scales up to a value of about 3.0 at larger values of a, and then starts to fall for a 

greater than about 400. Fitting the data to a function of the form D(a) = o H we obtain 

c + a 

/ X 10.51 , , 

D(c7 =3.03-— — (26 

^ ^ 17.87 + a ^ ^ 

Fitting the data to D(a) = a + he~^'^ , as suggested in [TT], we obtain 

D{a) = 3.19 - 0.97e-°-°^3'' (27) 

From this, we get -D(O) = 2.22 and D{oo) = 3.19, which agrees well with -D(O) = 2.12 and 
D{oo) = 2.98 reported in [Jl]. 

The smooth reduction of the dimension of the spacetime at shorter scales is one of the most 
interesting features of the CDT model. This is an indication of the highly non-classical behavior at 
small scales. The reduction in dimension for larger a is due to the finite (and rather small) volume 
of the spacetimes being considered. 

The spectral dimension vs step size plot for a (3+l)-dimensional spacetime is shown in figure |9j 
The spectral dimension in this case ranges from a value of about 1.8 at short scales up to a value 
of about 4.0 at larger values of cr, and then starts to fall for a greater than about 600, with the 
functional form being given by: 

This is to be compared with the 



reported in [2]. 



4.3 Spectral Dimension of the Spatial Slices 

The spectral dimension of a spatial slice can be computed using a process identical to that used for 
determining the spectral dimension of the full spacetime. We select the slice with the maximum 
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Figure 9: Spectral dimension for a (3+1) spacetime vs cr, = 80000, kq = 3.6, A = 0.6 



(d — l)-volume and initiate diffusion processes from a randomly selected (d — l)-simplex on this 
slice. 

The spectral dimension of the maximum volume slice for a (3+l)-dimensional spacetime is show 
in figure [Toj For short scales, up to o" ~ 40, the return probability differs considerably for even 
and odd step sizes, and the plot bifurcates. Once the even and odd step size spectral dimensions 
converge, D{a) is approximately 1.5, before a gradual reduction due to finite size scaling. In [1], 
a spatial slice spectral dimension of 1.56 is reported. This is value is computed by calculating the 
average return probability of all the spatial slices. The agreement of our results with value reported 
in [Ij indicates that each spatial slice has roughly the same spectral dimension. 



Figure 11 shows a corresponding plot for (2+l)-dimensional spacetime. In this case, the even 

150. The spectral dimension value then 



and odd a plots converge at a much larger value of a 
levels off at about 1.65 



Volume Profile 

At the beginning of this section, we discussed the emergence of an extended phase, which we said 
was a solution to the classical Einstein equations. However, the only evidence we have presented 
so far are figures [3] and [7j where the shape of the volume profile is indicative of such a solution. 
There is a rigorous, quantitative procedure, in which volume profiles of the spatial slices are fit to 
the Euclidean de Sitter solution, in 2+1 and 3+1 dimensions. This procedure is described in detail 
in [U]. We have been able to verify these results in 2+1 and 3+1 dimensions; our analysis will be 
presented in [5]. 



5 Conclusions 



The Causal Dynamical Triangulations approach to quantum gravity has many promising features, 
and our computer simulation has been able to successfully reproduce these results in 2+1 and 3+1 
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dimensions. Most of these results were first reported by Ambj0rn, Jurkiewicz, and Loll [Bill OH], 
and our simulation provides the first completely independent verification. We observe the emergence 
of a well defined extended geometric phase. This geometric phase shows a smoothly changing 
spectral dimension that ranges from a value close to 2.0 to a value larger than 3.0 for the (2+1)- 
dimensional spacetimes and from about 2.0 to 4.0 for the (3+l)-dimensional spacetimes. The 
spatial slices exhibit a non-classical behavior, as indicated by their spectral dimension plots. 

Our successful reproduction of many of the results reported by Ambj0rn, Jurkiewicz, and Loll, 
using a completely independent software implementation, should improve confidence in the CDT 
model. With the exception of the value of the phase transition point in the (2+l)-dimensional 
model, our results are in good agreement with those reported in [3l HI [H], and convince us of the 
overall validity of the CDT approach. 

Having successful reproduced the basic results of the CDT model, we are now using our CDT 
implementation to explore some further interesting problems. These include investigating the quan- 
tum fluctuations in (2+l)-dimensional CDTs and comparing the results with canonical quantization 
approaches [6], detailed analysis of the volume profiles in 2+1 and 3+1 dimensions [5] and explor- 
ing Horava-Lifshitz gravity on the CDT lattice [_7j. In the near future, we plan to release our 
implementation into the public domain. 
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